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Abstract 

Two different physical systems are said to be thermodynamically equiv¬ 
alent if one of the thermodynamic potentials of the first system is pro¬ 
portional to the corresponding potential of the second system after 
expressing the state variables of the first system in terms of those of 
the second by a transformation reversible throughout the state pa¬ 
rameter domain. The thermodynamic equivalence has a transitive 
nature so that physical systems divide into classes of thermodynam¬ 
ically equivalent systems that have similar phase diagrams. A first 
class of thermodynamically equivalent systems is formed by the ideal 
classical and quantum Fermi gases, whatever the dimensions of the 
confining spaces, and the one dimensional hard rod gas. A second 
class is formed by the physical systems characterized by interactions 
that coincide by a scaling of the distance and the coupling constant. 
A third class is formed by the ideal Boses gases in arbitrary spatial 
dimensions. The thermodynamic equivalence can also be dehned in a 
more general way. By so doing the first and the third class combine 
into a single class of thermodynamically equivalent systems. 

Keywords: Statistical thermodynamics, thermodynamic equivalence; 
hard-rod fluid; ideal classical gases; ideal quantum gases. 
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1 Introduction 


The thermodynamic behaviour of an open physical system is fully determined 
by the knowledge of one of its thermodynamic potentials which depends on a 
triple of state variables (one of which, at least, extensive) chosen among the 
three pairs of conjugate variables (V, p), {S, T) and (iV, p). Here, the exten¬ 
sive variables V, S and N respectively denote the volume, the entropy and 
the particle number of the sample while the intensive variables p, T and p are 
the pressure, the temperature and the chemical potentialt^’^^. For instance, 
choosing V, S and N as state variables, the thermodynamic potential that 
fully determines the thermodynamic behavior of the system is the internal 
energy U(y, S, N) equal to 

U{V, S, N) = -p{V, S,N)V + T{V, S,N)S + p{V, S, N) N (1) 

with 


P = p{V,S, N) 
T = T{V, S, N) 
p = p{V, S, N) 


dU{V, S, N) 

an 

dU{V, S, N) 

as ’ 

dU{V, S, N) 

m ■ 


( 2 ) 

(3) 

(4) 


If one wishes to determine another thermodynamic potential, say the grand 
potential T, p) = —pV, one solves equations ([3]) and (jl]) with respect 
to S and N so as to get 


S = S{V, T, p) and N = N{V, T, p). (5) 

After substituting these into equation([2]) one hnds the expression of the grand 
potential in terms of its natural variables, i. e. 


T, p) = -p(H, T, p) V. (6) 

In this way the transformation (V, S, N) (H, T, p), dehned by 

V^V, S ^ S{V, T, p) and N N{V, T, p), (7) 


allows one to describe the thermodynamic behaviour of the system using the 
grand potential instead of the internal energy. 

Consider now two different physical systems A and B in thermodynamic 
equilibrium. One wonders on the implications of a reversible transformation 
between the state variables relevant to a thermodynamic potential of the 
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first system and the state variables of the same thermodynamic potential 
relevant to the second system. These transformations lead to the dehnition 
of the thermodynamic equivalence (TE) that in turns unihes the thermody¬ 
namic behaviour of different physical systems. To the author knowledge, 
Lee^^l hrst spoke of TE when he realized that the grand potentials of the 
ideal Fermi and Bose gases are related to the same mathematical functions. 
Moreover, after equating the temperatures and the volumes of the two gases, 
Leed-5] got, in the two dimensional case, the fugacity transformation that 
converts the grand potential of the Fermi gas into that of the Bose gas plus a 
known function. Some years later, Anghel^®] showed that the last result also 
applies to two dimensional systems with particles obeying fractional statis¬ 
tics. This paper exploits the consequences of the TE dehnition that, as shown 
in section 2, can be formulated into a restricted and a generalized way. The 
analysis will be mainly conhned to the hrst dehnition because this ensures 
that the state diagrams of thermodynamically equivalent {te) systems are 
similarly shaped. Three examples of diherent physical systems that are te 
are discussed in sections 3 (ideal classical gases), 4 (systems characterized 
by a scalable interaction) and 5 (quantum ideal gases). The last section also 
shows that two physical systems, which are not te if one adopts the restricted 
dehnition of TE, become te with the more general dehnition of TE. Finally, 
section 6 draws the hnal conclusions. 


2 Thermodynamic equivalence 

For dehniteness let us refer to the grand potentials r2yi(l/A, T^, /i^) and 
Tb-, hn) of systems A and B and denote the transformation from 
(14, Ts, /is) of B to (14, Ta, /iA) of A as 


Vb ^ Va = Tb, /is) = VBVa,biTB, /ifi), (8) 

Tb ^ Ta = T^^g{VB, Tb, /is) = ta,b{TB, /is), (9) 

t /XA = hs) = Jta,b{TB, /^b)- (10) 

It is noted that the functional forms reported in the rightmost members of 

the above three relations are dictated by the extensive or intensive nature of 
the considered variables. The transformation is required to be endowed of 
continuous partial derivatives at least up to the second order, to be reversible 
and to be physically meaningful. The last condition implies that Va^iTB, /in) 
and tafiiTB, Hb) must be strictly positive and the second that the codomains 
of functions ([8])-([9]) coincide with the physical ranges of variables 14, Ta 


3 


and Ha- Denoting the inverse functions of ta^biTs, fin) and fia,b{TB, /ie) as 
h,a{TAi /i^) and fib,a{TA, the inverse transformation of equations (|H])- 
ffTOj) reads 

Va ^ Vb = VB,Ai^A, Ta, ^a) = VA'Vb,a(TA, Ha), ( 11 ) 

Ta^Tb = Tgj^{VA, Ta, Ha) = tb,a{TA, Ha), (12) 

Ha t hb Hb,a^^a, '^a, Ha) — Hb,aiTA, Ha), (1^) 

where 

Vb,a(TA, Ha) — ^/Va,b(TB, Hb) ^ / Va,b(tb,aiTA, Ha) , Hb,a(TA, Ha))- (1"1) 

Applying the state variable transformations l|5 |) - ffTU]) to the grand potential 
of system A, one obtains the new function ^ dehned as 

^a,b{^b, Tb, Hb) = (15) 

^A{V^^g{VB, Tb, HB),T^g{yB, Tb, hb), Ha,b{^b, Tb, Hb))- 

Systems A and B are said to be thermodynamically equivalent {te) if function 
fi^B is proportional to the grand potential of system B, i.e. 

^A,Bi^B, Tb, Hb) = Cj^^VLBiVB, Tb, Hb), (16) 

being a real positive constant. 

From this dehnition immediately follows that the thermodynamic equiva¬ 
lence (TE) may occur if; a) the codomains of the two grand potentials VLa 
and Qb coincide by a scale transformation, and b) the state diagrams have a 
similar shape in the sense that the state diagram of system A converts into 
that of system B by the transformation (l8|)- ffT0|) . [In fact, the state diagram 
of A is determined by the state variable points where the second order partial 
derivatives of are discontinuous. Equality flT^ and the assumed continu¬ 
ity of the partial second order derivatives of the state variable transformation 
imply a discontinuous behaviour of the relevant second order derivatives of 
Db-] Besides, the TE dehnition implies that: c) if physical system A is 
te to B and this is te to system C, then systems A and C also are te. In 
fact, the TE of B and C implies the existence of a reversible transformation: 
Vc ^Vb = Vg ,., Tc-^ Tb = Tg^ and p,c ^ Hb = Hb,c [where Vg^, Tg^ 
and Hbc denote functions that depend on Vc, Tc and hc] such that 

^B{Vgc,Tgc,HB,c) = ^B.d^c, Tc, He) = Cg^VLc{Vc, Tq, Hc)- (17) 
On the other hand, the state variable transformation 



= V 

A,C 

q 

hi” 

0 

III 

Hb,c)^ 

( 18 ) 

Tc^Ta 

= T 

A,C 

iV T 

— A,B \ ^ B,C ’ B,C ’ 

Hb.c)"! 

( 19 ) 

He Ha 

= h^A.C 

q 

hi” 

0 

III 

h'B.c) 

( 20 ) 
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existes and is reversible. Applying this transformation to and using equa¬ 
tions flTBD and ffT7|) . one gets 


C^B^BiyE, Tb, fin) 


= ( 21 ) 


^B—Vg fJ ,TB—Tg u tfJ'B—fJ'B Q 


^B—Vg Q ’ ^B—Tg (^ , fJ-B—t^B U 

= ^AB^Bc^ciyc, Tc, /ic), 


which proves the TE of A and C if A is te to B and B to C. This property 
implies that the totality of the physical systems splits into different classes 
of thermodynamically equivalent systems. 

The TE also implies that each thermodynamic quantity of system A is re¬ 
lated to the corresponding quantity of system B in a known way. Just to 
give an example, consider the entropies Sa and Sb of systems A and B. 
They are given by the relations Sa = —dflAiVA, Ta, fiA)/dTA and Sb = 
—d^BiVB, Tb, ^B)/dTB- Taking the Tb derivative of equations ffT5|) - ffT6|) 
and recalling that the pressure and the system particle number are respec¬ 
tively given by 


p(E, r,/i) = -dn{V,T,fi)/dV, (22) 

N{V,T,fi) = -dQ{V,T,fi)/dfi, (23) 

one finds 

dt 

Sb(Vb, Tb, Bb) = Sa(Va, Tb, ba)^ + (24) 

PAiVA, Ta, BaWa ^ + ITa(Va, Ta, Ba)^, 

SO that the entropy of system B is a linear combination of the entropy, the 
pressure and the number particle of system A and the coefficients of the lin¬ 
ear combination are known since they are appropriate derivatives of the state 
variable transformation. The same property holds true for the pressure and 
the system particle number. Similar inter-relations exist for other quantities 
as the specihc heats at hxed particle number and at constant volume or at 
constant pressure, even though they become more involved since these quan¬ 
tities are related to higher order derivatives of the grand potentials. 

It is observed that the validity of equation ffTHj) is not sufficient to determine 
the coordinate transformation. In fact, if one tries to determine the transfor¬ 
mation by a series expansion around a particular thermodynamic state, one 
realizes that the number of the involved unknown derivatives is larger that 
the number of the equations. However if one identihes two of the state vari¬ 
ables by setting, say, Va = Vb = V and Ta = Tb = T, then equation flT6|) 
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might be sufficient to express in terms of /i^ and T by a series expansion 
an, in this way, to determine the state variable transformation that converts 
into VIb- It niust be noted, however, that the resulting transformation 
as yet does not imply that two physical systems are thermodynamic equiv¬ 
alent. For this to happen, one must check that the resulting transformation 
bi-injectively maps the state variable domain of A onto that of B. 

It is also noted that, in dehning the TE, one has to choose a thermodynamic 
potential. One wonders if, assuming that two systems are TE with respect 
to, say, the grand potentials, they also are TE with respect to another ther¬ 
modynamic potential. In section 3 it is reported a case where this happens, 
but the general answer is as yet unknown. 

Finally, the reported dehnition of TE could be made less restrictive modifying 
condition flTB]) as follows 



where and are known functions of the reported state variables. 
This new dehnition, in contrast to that given by Eq. ([22]), will be referred 
to as generalized thermodynamic equivalence (GTE). It obeys the transitive 
property [he. c)] but does not obey properties a) and h). Unless one restricts 
the analytic form of Aabi the GTE dehnition looks too much general to 
become nearly trivial. In fact, any invertible transformation of the state 
variables of systems A into those of B would yield the GTE of A and B 
by setting Aab = (Has — ^b) and Cab = 1- For this reason the following 
analysis will be mainly conhned to dehnition ffT6|) . 

3 The hard rod gas in one dimension and 
the ideal classical gas in any dimension 

This section reports a hrst set of physical systems that are thermodynami¬ 
cally equivalent. They are the classical ideal gases in any space dimension D 
and the hard rod gas in one dimension. Gonsider hrst the pair formed by the 
one dimensional and the three dimensional classical ideal gases. The grand 
potential of the classical ideal gas, conhned into a box of length side L and 
contained into a space of dimension D(= 1,2,3), reads^^’^l 



(26) 


with Vg = and 


\ = h/ \/2TTmkgT 


(27) 
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denoting the de Broglie length, h and kg the Planck and the Boltzmann con¬ 
stants and m the particle mass. For later reference it is convenient to report 
here also the expressions of the particle nnmber density and the pressnre, 
respectively given by 




1 ^ 

Vo AO 


(28) 


and 

Tn„(T, ft,). (29) 

To prove the TE of the considered ideal gases, one writes the right hand side 
of eqnation fl26l) nsing dehnition fl?7|) as 


kg{27rm)^/^ u,,,o/kgTg 

J^D D V gC. 


It is clear that the transformation 


(30) 


n = n/io. 


rri. 

hig, 3-^; 


. 2/3 


(31) 


where /q denotes an arbitrarily chosen length, is invertible and continnons 
with all its partial derivatives. Once it is applied to ^(Ei,Ti,/ij^ i) one 
hnds that 


k'ig, - 


2tt m kg 


12 ^ig, 3 (^ 3 ) ^ 3 ) k'ig,3)i 
'■0 


(32) 


which is eqnation (IT^ with C^^g = /27imkgll. Thns, the TE of the one¬ 
dimensional and the three-dimensional classical ideal gases is proved. By 
the same procednre one proves the TE of the one-dimensional and the two- 
dimensional ideal gases. In this way, by the transitiveness of the TE [see 
property (c) below eqnation ffT6|) ]. one conclndes that the D-dimensional 
ideal classical gases, whatever the positive integer D valne, form a class of te 
systems. 

To prove that the one-dimensional ideal gas is te to the one-dimensional 
hard rod gas of particles of mass m and length a, it is hrst recalled that 
the last system is characterized by the interaction potential (r) dehned as 
Vhr{f') = 0 if r > cr and v^^{r) = cx) if 0 < r < a. Rayleighand Tonk^®] have 
since long determined the eqnation of state of the one dimensional hard rod 
system. More recently Robledo and Rowlinsonl®] have determined the fnll 
thermodynamic behavionr of this system within the grand canonical frame¬ 
work paying great attention to the behavionr of the one particle distribntion 


7 







function close to the system walls. The mentioned TE between the ideal gas 
and the hard rod gas holds only true in the thermodynamic limit, i.e. when 
the walls are inhnitely far away. The grand partition function of the hard 
rod gas, conhned to a box of size L, (after correcting a misprint present in 
the corresponding expression of reference [9] ) is 


A k 


IhriL, T, ^ ^ 


k =0 


(33) 


with 


C,, = (34) 

and [L/a] equal to the integer part oi L/a. The mean particle number of the 
particles (N) is 


(Af), = c T - tr - tr (iV>,_ J. 


SO that, in the limit L — )■ cxo with (N)^/L = Uhr fixed [uhr denoting the bulk 
particle number density of the hard rod fluid], one hnds that 


lim C 


CT, T,/i^^) 


L—)-oo 


"-^hriL, T, 


1 - an. 


It also results that 


5 log! 


= c. 


-‘hr 


(L-CT,r,fi,,) 


dL Si,(L,r,ft„) 

Using equation fl36|) the above relation can be integrated to yield 

+ const, 

-L (J IL, 

fir 

which in turns implies that 

Ehr{L - a, T, ^J/Ehr{L, T, = exp(-an^7(l - anj). 
This result converts equation fl36l) into 

n 


(36) 


(37) 


(38) 


Chr 


1 - an. 


hr ^ 1 — n, a 
6 


(39) 


that allows one to relate the chemical potential to the particle number 
density n^^. Putting const = 0 in equation ([3H]) to ensure the extensiveness 
of the function, one hnds that the grand potential of the hard rod system is 


= -/c^Tlog(Sfc^(L,T,p^J) = -k^T- 


Ln, 


crn^ 


(40) 











This expression depends on while it is important to know in 
terms of in order to derive from it the other thermodynamic quantities. 
This can easily be realized using equation fl3^ to express in terms of . 
To this aim one recalls that the solution, with respect to x, of the equation 
y = xe^ is a transcendental function known as Lambert’s functiont^^’^^i and 
denoted as x = W{y). This function also is the solution of the following 
differential equation 


W\y) 


W{y) 

y{l + W{y)) 


(,-W{y) 

1 + W{y) 


(41) 


with the boundary condition 1T(0) = 0. Hence, the solution of equation fl3^ 
with respect to is^^^l 




Wiae^hrl^B'^ jX) 
cr(l + Wiae^hrl^a'^ jX)) 


(42) 


Recalling that W{y) ~ log(|/) as ?/ —)■ oo, this equation implies that —?■ 
1/(T as —)■ oo and 1/a clearly represents the system density at the closest 
packing. The substitution of fl42|) within equation fldH]) yields the expression 
of the hard rod grand potential in terms of its natural variables, i. e. 

njL,T,i,,J = ( 43 ) 

[The use of the same symbol here and in equation (ITU]) is dictated by the 
sake of notational simplicity even though it is not mathematically correct. 
This same convention will be used in a few of other points later.] The inter¬ 
action of the hard rod fluid vanishes in the limit a —)■ 0 so that the hard rod 
fluid must behave as the ideal gas in the aforesaid limit. Since the Lambert 
function is such that W{y) ^ y as y ^ 0 one immediately realizes that 

= -k.TLe^^rV/x = ( 44 ) 


once one sets and the identical behaviour of the two systems in 

the limit a —?■ 0 is proved. 

The interesting point is that the hard rod and the ideal gases are fe if a > 0. 
Consider, in fact, the following transformation of the state variables of the 
two systems 


T=T=T, L. =L=L (45) 

zg hr ' zg hr \ / 

= kgT{\og{\/a) + \og{W{\e^^hrl^BT 
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It is a reversible transformation because the last equality can be inverted 
throughout the full range (—cxo, cx)) of to yield 

k T 

+ (46) 

Once transformation (H5|l is applied to equation (l26ll one hnds that 


k'igi'^1 k'hr)) ~ ^hri^l'^lk‘hr)l 


(47) 


that is equation ([TB]) with = 1. In this way the TE of the one dimensional 
hard rod and ideal gases is proved. Due to property (c), one concludes that 
the D-dimensional ideal gases and the one dimensional hard rod gas form a 
class of te systems. 

The substitution of fITS]) into equation fl25D and the use of equation flT2D imply 
that the particle number densities of the two systems are related by 


n 


■^9 


n 


hr 


1 — an, 

hr 


( 48 ) 


This equation makes it evident that, as n^^ ranges over the allowed physical 
range [0, 1/cr], n. ranges over its physical range, namely the positive half¬ 
axis. The equation of state of the hard rod system is obtained by taking the 
derivative of with respect to (—T). It can directly be expressed in terms 
of variables T and derivating equation fldH]) . One hnds 


Pnr 


, ^ Ln. 
k T -^ 


cm. 


( 49 ) 


and, once one uses here relation (1481) . it coincides with that of the one¬ 
dimensional ideal gas p. = k^Tn. . 

The entropy expressions are easily obtained from the grand potential expres¬ 
sions. For the hard rods huid one hnds 


S,,{L,T,gJ 


^ir(ffe'v/<*i>T)/A) 

a 


-3 

.2 


+ W{aef^^r/ikBT)/x) 


hhr _i_■ 

kgT 1 + W{ae^hr/('^B'^'i / Xy’’ 


which, in terms of reads 


S„,{L,T,n,_^) = i, Ln,,(3/2- log(AnjJ(l -n,,<T))). 


( 50 ) 


( 51 ) 


For the one-dimensional ideal gas one gets 


S.JL,T,g,) 


kg Lm\/n/ b'^(^ kgT12 — 2p.^/{h\/kgmT^ 
^B^’^,s(3/2-log(An,j). (52) 
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The substitution of transformation fHS]) into flS2]) does not convert the ideal 
gas entropy in that of the hard rods. In fact, the two entropies are related 
by the relation [see equation fl24|l ] 


because the /i.^ transformation, given by fl45hh depends also on T. 

The expressions of the Helmotz free energies and internal energies of the two 


systems are now reported for completeness. The free energies are 



= -k^ 

Tnr (l - log(An^7(l - n^^a))^ 

(54) 

and 


= -*:,r„iV„(l-log(Anj). 

(55) 

and the internal energies read 




UJL,„S,„NJ = 

^ hhr ^hr ^ — 3+2Sf^^/(kpN^^) 

(56) 


Arm (1 — 

and 

U. (L. ,S. ,N. 

tg \ tgi tgi tg 

Anm 

(57) 


One easily verihes that the following reversible transformation of the state 
variables 


=K =N, T. =T =T, L = L. (L, ,N) = L, -Na (58) 

ig hr ' tg hr ’ %g *9 V hr ' / hr \ / 

(with the bound N < L^_^/a) yields 

f’.,(r„(U..A'),T,Af) =F,,(L,„T,Af). (59) 

This relation proves the TE of the two systems making use of the Helmotz 
free energies. 

If one considers, as thermodynamic potentials, the internal energies [see equa¬ 
tions 0561) and 0571) ]. it is straightforward to show that the following reversible 
state variable transformation 

77 

L. = L = L, n. = - — —-;r 7 ^, S. In. = S jn 

1-9 hr 1 tg ^ _ (JTl i-9 > i-9 hr / hr 

yields 

= U,AK,S,^,NJ, (60) 
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i.e. the internal energies also are te. Eqnations and fl^ snggest that if 
the grand potentials are te the same happens to the the other thermodynamic 
potentials, even thongh a general proof of this conjectnre is still lacking. 
Before conclnding the section, we report two properties of the hard rod flnid. 
The hrst concerns the heat capacities at hxed particle nnmber and at hxed 
volnme or at hxed pressnre. The hrst is easily derived from expression (15 ip 
and the second from the same expression after expressing in terms of 
by eqnation (H^ . In this way one respectively hnds 

^nd C (61) 

which have the same analytic forms of the ideal gas heat capacities, i.e. 
C. k„N. /2 and C. = 3k„N. /2. 

The second property concerns the existence of the proportionality relation 
between the grand potential and the internal energy, i.e. 

S,„ NJ = NJ (62) 

-L (J IL, 

nr 

This relation is obtained by snbstitnting the temperatnre expression T = 
—dU/dS, resnlting from eqnation (15Bp . into eqnation (HUp and then by com¬ 
paring the resnlt with ([56]). For the ideal gas the eqnivalent relation is^^^^ is 
ff.g j that dihers from (l62p because the proportionality coefficient 

is here strictly constant while it depends on the particle number density for 
the hard rod gas. 


4 Systems having the same scaled interaction 


The systems considered in the previous section are characterized by interac¬ 
tions that coincide if the inter-particle distance exceeds a. In this section one 
considers the cases where the interactions have a similar behaviour through¬ 
out the full range of distances in so far they have the form gV{r/a) with V (r) 
hxed and g and a variable. Hence, one assumes that system A is made up of 
particles of mass m^, diameter a a, and interaction gAV (j" / (^a) and system 
B of particles of mass nig, diameter as and interaction gsVir/aB)- The 
two systems are te. In fact, the classical grand canonical partition function 
of system A, after integrating over the particles’ momentums, is 


i{VA,TA, Pa) = ^ 


za^ {2'Kni^kBTAy 


A:=0 


klh^^ 


X 



JA 




'Va 


(63) 


12 





with 




The corresponding grand potential is 

QA{VA,TA,flA) = -kBTA\og{E^{VA,TA,HA))- 
The grand partition function of system B is 


^-‘siyByTB, ^ 


zb’‘ {27imgkBTBY 


k=0 


k\h^^ 


X 






’Vg 


'Vg 


(64) 

(65) 


( 66 ) 

(67) 


Putting aAB = (^aI'^b (and = Ogjo^ and performing the following 
change of the integration variables r* —)■ a ab'^' ji conhguration integral 

present in the /cth term of series becomes 


(y AB 


3k 


' kgTA Si<i<j<fc n(r-'ij7(Tg)^7Vj,/ 


' VAO-f 


' VAO-f 


At this point one easily verihes that the state variable transformation 
Vb = Vb{Va) = VACJg^Y Tb = Tb{Ta) = Qb Ta/qa, 

and Zb = za- 

the last being equivalent to 




2 

fiB = g^BiTA, Ha) = — yUA + -fcijTA log 
9a ^ ^ 


^aQa^A 
mg gg(Tg^ 


is reversible and yields 

'^B{VB(yA),TB(TA), Hb(Ta, Ha)) = '^A(yA,TA, Ha)- 
This, in turns, implies that 

klB{VBiVA),TB{TA), Hb{Ta, Ha)) = igg/gA)^A{VA,TA, Ha)- 


( 68 ) 


(69) 


(70) 

(71) 

(72) 


In this way, according to equation flTBD with Cg^ = gg/gj^-, the TE of systems 
A and B is proved. This result shows that a class of equivalent thermody¬ 
namic systems is formed by those systems such that each of them is formed 
by particles, of a given mass and diameter, which interact with a potential 
having a particular coupling value and a hxed shape. Hence, the systems 
interacting with a Lennard-Jones potential form a class of te systems. The 
same happens for the Coulombian systems. It is noted that the systems, 
which are te under the variable transformations fl69ll and (1701) . also obey the 
law of corresponding states^^^. 
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5 The ideal quantum gases 


It is now shown that, within a space of dimension D, the ideal classical gas 
is te to the ideal Fermi one and is not te to the ideal Bose one if one adopts 
equation ffT^ as the TE dehnition. To this aim it is hrst recalled that Leet^^ 
showed that the grand potential expressions of the two ideal quantum gases 
can be expressed in terms of the same mathematical functions as follows 


^B,d{^B,Tb, fig) 
^F,Di^F,TF, 


9b VshgTi 


9f kgTp 


Li_D/2+l(^_B), 

Liz)/2+i(— 


(73) 

(74) 


Here suffices B and F respectively refer to the Bose and the Fermi gas, Qg = 
2sg +1 with Sg equal to the particle spin, X = h/^2'KmkgT, Vb = Vp = 
with L equal to the length of the conhning box edge and D equal to the 
space dimensionality. Further, zb{= and zf{= are the 

Bose and Fermi gas fugacities. Finally, Lis(z) is the so-called polylogarith- 
mic functiont^^] of index s that is simply relatedl^"^] to <F(z,s,a), the Lerch 
function, by the relation 


Lh(z) = z<F(z,s, 1), 


(75) 


so that the ideal Fermi and Bose grand potentials can also be expressed^^^^ 
in terms of the Lerch function. This is dehned as^^^l 


$( 2 ;, s, a) = 


00 ^5 — 1 ^ — 


F(5) 


z e 


—t 


dt 


(76) 


for Re(a) > 0, Re(s) > 0 and Re(. 2 ) < 1. It obeys the three properties 




(P + “)■ 


if l^l < 1, 


(77) 


d 

^{z,s-l,a) = {a +z-^^{z,s,a) (78) 

and 

1 d 

<h(z, s1, a) = ——$( 2 ;, s,a). (79) 

s oa 

The integral expression 0761) shows that the Lerch function is analytic in the 
complex z plane cut along the real axis from one to 00 . Further, equation 075]) 
shows that $( 2 ;, s, a) > 0 as 2 ; ranges over (—cx), 1). From this property and 
equation 075|) it follows that l-X{z) >0if0< 2 ; <1 and Iji{z) < 0 if 2 ; < 0. 


14 












Recalling that Li'( 2 ;), the derivative of Li( 2 ;), obeys Li'( 2 ;) = ljis_i{z)/z [this 
relation immediately follows from equations fl7S]) and fl75]) ]. it follows that 
Li; {z) > 0 throughout (— 00 , 1). One concludes that Lis(z) (with s > 1) 
monotonously increases as z goes from —00 to 1. The behaviour of the 
LL(z)s for s = |, 2, |, the values associated to D = 1, 2 and 3, are shown 
in Figure 1 for — 00 < z < 1. Further, the leading behaviour of Lis(z), as 
z -)■ 1“, ist^^l 

Lh(2;) ^ C(s) - C(s - 1)<5+^[C(s-2) - C(s-+ (80) 

(5"-^F(l-s)-^(s-l)(5T(l-s), as 5 = (1 - z) ^ 0+, 

where ({x) is the Riemann zeta function. As ^ cxo, Lis(— 2 ;) behaves as^^’^l 

log* \z\ 

Lh(- 2 ;) Ri - as z ^ cx). (81) 

F(1 + s) 

[This behaviour is obtained using Joncquiere’s relation (see equation (1.11.1.16) 
of reference [14]) and the asymptotic expansion of the Hurwitz function C,{a, b) 
as 6 ^ cxo in the complex plane (see equation (25.11.43) of reference [10])]. 

At this point, for each D value, the proof of the TE of the classical ideal 
gas and the ideal Fermi gas is straightforward. In fact, consider the state 
variable transformation 


T.^=Tf = T, Vig = VF = V and Zig = -Ud/2+i{-zf), (82) 

the last relation being equivalent to 

ft, = ft,(AF,T) = A^.riog(-LiD/ 2 +i(-eWft.^)). (83) 


It is dehned throughout the physical ranges of variables Tp, Vf and zf and the 
ranges of the resulting variables T.^, Vig and /r.g coincide with their relevant 
physical ranges. The transformations is also reversible. [This property, trivial 
for the hrst two variables, holds true for variables Zig and zf owing to the 
reported properties of the relevant Lis(z) and it is also made evident by 
Figure 1.] Then, the substitution of fl52D into fl26]) and the comparison of the 
result with equation fl74|) yield 




( 84 ) 


that is condition flTHl) with = 1/gp. A corollary of this result is the 
property that, whatever the positive integer D, the ideal Fermi gases and 
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Li,(z) 



Figure 1: The curves in colours are the plots of Li 5 ( 2 ;) for s = 5/2 (blue), 
s = 2 (red) and s = 3/2 (magenta) within the range —oo < z < 1. The 
curves are broken within the negative ^ half-axis and continuous within the 
range 0 < z < 1 to make it more evident that they respectively are the plots 
of Lh(—Zi?) with zp > 0 and of IAs^zb) with 0 < < 1. Note that the 

Li 5 (z) limits are hnite as z ^ 1. The black broken and continuous straight 
lines are the plots of —Zig and Zig as /i.^ runs over (—cxo, cx)). The lower 
horizontal thin line shows that the equation —Zig = Lis(— 2 ;i?) has a single 
root for each Zig or zp value as well as for each of the considered s values. 
On the contrary, the upper horizontal thin line illustrates a case where the 
equation Zig = IAs^zb) has no roots. 

the ideal classical gases as well as the one dimensional hard rod gas are all te 
among themselves, i. e. they form a class of te systems. This property follows 
from the transitiveness property (c) and the TE of the D-dimensional ideal 
classical gases proved in section 3. 

This result might look surprising because the classical ideal gas, in contrast 
with the Fermi one, is physically inconsistent over a much wider range of the 
state variables as it happens, e.g., for the entropies and the specihc heats. 
Nonetheless no contradiction is possible. In fact, it was already noted that the 
TE does not imply the equality of corresponding thermodynamic quantities. 
Consider, for instance, the entropy case. The entropy of the Fermi gas is 
always positive^^^ while that of the ideal gas can be negative!^®!. Since they 
are related as in equation i.e. 

Se(V, T, fip) = S„{T, V, 

it follows that the right hand side turns out to be everywhere positive thanks 
to the contribution related to the partial T-derivative of ^^^{T,fip). 

Finally, the ideal quantum Bose gases, compared to the ideal classical and 
quantum Fermi gases, are now discussed from the point of view of the TE. 
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Interestingly, the conclusion depends on the adopted dehnition of TE. Con¬ 
sider Erst the restricted dehnition of TE given by Eq. flTB]) and the ideal Bose 
gases with D = 1 and D = 3. Assume also that the temperatures and the 
volumes are related as in Eq. Then the expression within the square 

brackets in Eq. ([73]), in the case D = 1, becomes 


9b ^BikgTsi 

Ai 


-)■ 


9b ^bs^bTes 



I 


2 

0 


(85) 


that, aside for the factor coincides with the expression relevant to the 
case D = 3. The TE is ensured if one proves the existence of an invertible 
transformation zbi = zbi{zb 3 ) with 0 < zbi < 1 as zb 3 ranges over [ 0 , 1 ) 
such that 


^k/2{zBl{zB3)) 


Li3/2(1) 

Li5/2(1) 


^h/2{zB^ ■ 


( 86 ) 


Owing to the reported properties of the polylogarithmic functions, the above 



Figure 2: The blue and the red curves respectively are the plots of 113 / 2 ( 2 :) 
and Li 3 / 2 (l)Li 5 / 2 (.s)/Li 5 / 2 (l). Each horizontal line, of height smaller than 
113 / 2 ( 1 ) intersects the two curves at A and B. The abscissa of A is uniquely 
associated to that of B so as to dehne zbi = zbi{zb 3 )- 

equation uniquely determines zbi as function of zb 3 - Besides, the resulting 
function is continuous, endowed of continuos derivatives of any order, spans 
the interval [0, 1) and is invertible within this interval. Figure 2 makes these 
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properties evident. Substituting Eqs. fISS]) and flSB]) within fl75]) one obtains 




^^Li5/2(1) 


/o27rm/Cg, Li3/2(1)' 


^s,3 (^^3) ^3) hs.s) 


( 87 ) 


where ^ j. Comparing fl87)) to ffTHl) . one 

concludes that the ideal Bose gases in one and three dimensions are te. The 
proof of the TE of the ideal Bose gases for the cases D = 2 and D = 3 
proceeds quite similarly. Therefore, the ideal Bose gases, whatever the space 
dimension D, form a class of thermodynamic equivalence. 

This class is different from that of the ideal Fermi gases and the classical ideal 
ones. To prove this statement it is sufficient to show that one physical system 
belonging to the hrst class is not TE to a system belonging to the second 
class. Consider then the ideal classical gas and the Bose quantum one in the 
case D = 3 and assume that the temperatures and the volumes of the two 
systems are equal. One observes that, as zb ranges over the physical range 
[0,1), function ^ 5 / 2 ( 25 ) ranges over [0,(C(5/2)) [see equation (180)) ]. Hence, 


at hxed T, 


9b ^ 

\o 


ranges between [0,C(s)) as /i^ spans 


its physical range (— 00 , 0). On the contrary, for the ideal gas, one hnds that 

ranges over (—cxo, 0) as /r. ranges over its physical 






range {— 00 , 00 ). The ranges of the two grand potentials are different and, 
therefore, the necessary condition [see property a) of section 2] for the two 
systems to be te, in the restricted sense of Eq. (1T6|) . is not obeyed. Thus, 
the proof of the non equivalence is achieved. 

This conclusion completely changes if one adopts the GTE dehnition given by 
Eq. (123]) . In fact, it will now be shown that, if one adopts the GTE dehnition, 
the aforesaid two distinct classes of restricted te systems combines into a 
single class of generalized te systems. To this aim, it is hrst observed that 
the GTE implies the restricted TE and that transitive property c) of section 
2 applies to both dehnitions. Then one recalls the result obtained by Lee^^^, 
namely: the one-to-one fugacity transformation zp = zf{zb) = ^n/(l — zb), 
in the case D = 2, relates the grand potentials and the one particle densities 
of the Bose and Fermi gases (outside the Bose condensation region so as 
to avoid Pathria’s criticism^^^ and under the simplifying assumption that 
Qg = Qp = 1) a.s follows 






k^TV 


\\T)pg\zB), 


( 88 ) 


and 


PbP^^b) — Pp,2{zf{zb)) — ^Lii(zs). 


( 89 ) 
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The comparison of Eq. with fISSD and flS^ shows that the two-dimensional 
Fermi gas is GTE to the two-dimensional Bose gas with ^^^(1/,T,/i^) = 
kgTV\Ai^{zB)/y?■ In this way the GTE of the classical ideal gases as well 
as the Fermi and the Bose quantum gases, whatever the space dimension D, is 
fully proved. One should also note that the added contribution kgTVLii‘^{zB )/ 
ranges ove [0, oo) as zb ranges over [0, 1) so that the codomain of the right 
hand side becomes inhnitely large as it happens to the function reported on 
the left hand side. 


6 Conclusion 

The consideration of the reversible transformations between the state vari¬ 
ables of different physical systems naturally leads to the dehnition of ther¬ 
modynamically equivalent systems. This paper mainly considered the most 
restrictive dehnition of the TE which implies appropriate constraints on the 
codomains of the equivalent thermodynamic potentials as well as the sim¬ 
ilarity of the phase diagrams for the TE to occur. The most interesting 
feature of the dehnition is the possibility of dividing all the physical systems 
into classes of TE. A hrst class of te systems is formed by the ideal classical 
and quantum Fermi gases in any space dimensions and the one dimensional 
hard rod gas. A second class is formed by the systems where the particle 
interaction can be brought to coincide by a scaling of the distance as well as 
of the coupling constant. Since these systems show phase transitions, they 
are therefore endowed of similar phase diagrams. A third class is formed 
by the ideal Bose gases in any spatial dimension and this class is diherent 
from the hrst one. But if one adopts the generalized dehnition of the TE the 
hrst and the third class become a single class of generalized TE. It is also 
stressed that the TE of two diherent physical systems does not ensure that all 
the thermodynamic quantities of one system are equal to the corresponding 
quantities of the other system. It does only ensure that the quantities are 
linearly related among themselves with coefficients that depend on suitable 
(partial) derivatives of the state variable transformation. 
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